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Preface 


This  report  is  a  user's  guide  for  computer  program  CSMT  which  can  be 
used  to  calculate  shear,  moment,  and  thrust  for  a  given  structure  section  from 
stress  results  of  a  two-dimensional  finite  element  analysis.  The  work  in 
writing  the  program  and  this  user's  guide  was  accomplished  with  funds  provided 
to  the  Automatic  Data  Processing  (ADP)  Center,  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  Vicksburg,  Miss.,  by  the  Engineering  and  Construc¬ 
tion  Directorate  (ECD)  of  the  Office,  Chief  of  Engineers,  U.  S.  Army  (OCE) , 
under  the  Computer-Aided  Structural  Engineering  (CASE)  Project. 

Specifications  for  the  program  were  prepared  by  members  of  the  CASE  Task 
Group  on  the  Finite  Element  Method  of  Analysis.  Members  of  the  task  group 
during  the  period  of  development  of  the  program  were  as  follows: 

Mr.  Tom  McGee,  Nashville  District  (Chairman  until  February  1983) 

Mr.  Dave  Raisanen,  North  Pacific  Division  (Current  Chairman) 

Mr.  Rich  Flauaus,  St.  Louis  District 

Mr.  Dick  Huff,  Kansas  City  District 

Mr.  Paul  LaHoud,  Huntsville  Division 

Mr.  Ed  Ailing,  U.  S.  Soil  Conservation  Service 

Mr.  Jerry  Foster,  Federal  Energy  Regulatory  Commission 

Mr.  Paul  Noyes,  Seattle  District  (joined  in  May  1983) 

Mr.  Lucian  Guthrie,  OCE 
Mr.  Robert  L.  Hall,  WES 
Dr.  N.  Radhakrishnan,  WES 

The  program  and  this  report  were  written  by  Mr.  Fred  T.  Tracy,  Chief, 
Research  and  Development  Software  (RADS)  Group,  ADP  Center,  WES;  Mr.  Robert  L. 
Hall,  Computer-Aided  Design  (CAD)  Group,  ADP  Center;  and  Mr.  Kenneth  W.  Trahan, 
summer  student  with  the  RADS  Group.  Initial  guild  line  for  input  data  was  pro¬ 
vided  by  Mr.  H.  Wayne  Jones,  CAD  Group.  The  work  was  managed  and  coordinated 
by  Dr.  Radhakrishnan,  Special  Technical  Assistant,  ADP  Center,  and  CASE  Project 
Manager.  Mr.  Lucian  Guthrie,  Structures  Branch,  ECD,  was  the  OCE  point  of 
contact . 

Commander  and  Director  of  WES  during  development  of  the  program  and 
preparation  and  publication  of  this  report  was  COL  Tilford  C.  Creel,  CE. 
Technical  Director  was  Mr.  F.  R.  Brown. 
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Conversion  Factors,  Non-SI  to  SI  (Metric) 
Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply _ 

feet 

kip  (1000  lb  mass)-feet 
pound  (mass) -feet 
pounds  (force) 

pounds  (force)  per  square  inch 


By _ 

0.3048 

138.255 

0.138255 

4.448222 

6.894757 


To  Obtain 
meters 

kilogram-meters 

kilogram-meters 

newtons 

kilopascals 
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USER'S  GUIDE:  COMPUTER  PROGRAM  TO 
CALCULATE  SHEAR,  MOMENT,  AND  THRUST  (CSMT) 

FROM  STRESS  RESULTS  OF  A  TWO-DIMENSIONAL 
FINITE  ELEMENT  ANALYSIS 

Purpose 

1.  The  purpose  of  program  CSMT*  is  to  calculate  shear,  moment,  and 
thrust,  for  sections  of  a  structure  specified  by  the  user,  from  stress  results 
of  a  two-dimensional  finite  element  (FE)  analysis. 

Program  Details 

Input 

2.  The  bulk  of  the  input  is  for  geometry  definition  (node  and  element 
data)  and  the  stress  results  from  a  FE  analysis.  The  node  and  element  data 
are  read  free-field  from  one  data  file,  and  the  stresses  are  read  from  another 
file.  The  remaining  data  are  interactive  commands  to  specify  section  informa¬ 
tion  and  to  obtain  grid,  section,  and  output  plots. 

Geometry  data 

3.  The  geometry  data  are  read  free-field  from  a  line-numbered  data 
file.  The  data  consist  of  the  node  and  element  data  which  define  the  FE  grid. 
Thus,  the  (X,Y)  coordinates  for  each  node  and  the  node  numbers  for  each  ele¬ 
ment  are  required.  Figure  1  shows  the  exact  data  file  format.  Note  that 
elements  may  be  4-  or  8-node  quadrilaterals. 

Triangular  elements 

4.  Triangular  elements  are  input  by  collapsing  one  side  of  a  quadri¬ 
lateral  to  a  point.  Figure  2  shows  3-  and  6-node  triangular  elements  and  the 
corresponding  element  data.  It  should  be  noted  that  all  input  is  the  same  for 
either  quadrilateral  or  triangular  elements. 


*  CSMT  is  designated  X0063  in  the  Conversationally  Oriented  Real-Time 

Program-Generating  System  (CORPS)  library.  Three  sheets  entitled  "Program 
Information"  have  been  hand-inserted  inside  the  front  cover  of  this  report. 
They  present  general  information  on  the  program  and  describe  how  it  can  be 
accessed.  If  procedures  used  to  access  this  and  other  CORPS  library 
programs  should  change,  recipients  of  this  report  will  be  furnished  a  re¬ 
vised  version  of  the  "Program  Information." 
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LN  NEL  I(NEL)  J(NEL)  K(NEL)  L(NEL)  Il(NEL)  Jl(NEL)  Kl(NEL)  LI (NEL)  M(NEL) 

LN  =  line  number 

NNP  =  total  number  of  nodes  in  grid 

NEL  =  total  number  of  elements  in  grid 

NNELE  =  number  of  nodes  per  element  (4  or  8) 

NMATP  =  number  of  material  properties 
X(N)  =  X  coordinates  of  node  N 

Y(N)  =  Y  coordinates  oi  node  N 

M(N)  =  material  type  for  element  N 


=  corner  nodes  of  an  element 


I1(N) 

J1(N) 

K1(N) 

L1(N) 


=  midside  element  nodes  (necessary  only  if  NNELE  =  8) 


Note:  Node  and  element  data  must  be  in  ascending  order. 

Figure  1.  Geometry  data  file  format 


tmiKaGn C 


ELEMENT 

DATA 


Element  Data: 

NEL  I  J  K  L  MAT 

54  10  21  35  35  1 

Material  Type  =  1 


a.  3-node  triangular 
element 


35 


b.  6-node  triangular 
element 


Figure  2.  Triangular  elements 


Global  stress  data 


5.  The  global  stress  data  are  also  read  free-field  from  a  line-numbered 

data  file.  These  data  consist  of  the  o  ,  a  ,  and  x  stresses.  These 

x  y  xy 

stresses  may  be  available  at  different  locations  in  an  element  depending  upon 
the  FE  program  used  and  the  type  of  element.  To  accommodate  the  various  pro¬ 
grams,  program  CSMT  accepts  the  following  options  for  computing  stresses: 

a.  At  the  "center"  of  each  element--Type  =  0. 

b.  At  each  node  (an  average  value  for  all  elements  meeting  at  a 
node) — Type  =  1. 

c.  At  each  node  for  every  element  (as  many  values  at  a  node  as  the 
number  of  elements  meeting  at  the  node) — Type  =  2. 

d.  At  each  midside  and  "center"  for  each  element  (for  4-node  ele¬ 
ments  only) — Type  =  3. 

Figure  3  shows  the  exact  data  file  format.  The  "center"  of  an  element  is 
assumed  to  be  at  the  coordinates  s  =  0  ,  t  =  0  for  a  parametric  coordinate 
system. 

Interactive  Commands 


6.  Interactive  commands  are  provided  to  give  the  user  flexibility  in 
defining  those  sections  along  which  moments,  shears,  and  thrusts  need  to  be 
calculated  and  in  obtaining  output.  These  commands  are: 

a.  SECTION.  Define  a  section. 

b.  PLOT.  Plot  the  grid  and/or  section. 

c.  WINDOW.  Do  a  window  plot. 

d.  OUTPUT.  Display  output  for  a  section. 

e.  END.  End  running  the  program. 

Because  all  the  commands  start  with  a  different  letter,  only  the  first  letter 
of  a  command  need  be  input.  However,  the  user  may  use  as  many  of  the  letters 
as  desired. 

SECTION 

7.  The  general  form  for  this  command  is: 

SECTION  I,  XI,  Yl,  X2,  Y2,  M 

This  command  allows  the  user  to  define  the  Ith  section  where  shear,  moment, 
and  thrust  calculations  are  desired  by  defining  two  points,  (XI,  Yl)  and 
(X2,  Y2),  representing  the  beginning  and  ending  points  of  the  section.  The  M 


LN  TYPE  TC  Control  Data 

STRESSES  FOR  TYPE  =  0 

LN  SX(1)  SY(1)  SXY(l) 


LN  SX(NEL)  SY(NEL)  SXY(NEL) 

STRESSES  FOR  TYPE  =  1 

LN  SX(1)  SY(1)  SXY(l) 


LN  SX(NNP)  SY(NNP)  SXY(NNP) 
STRESSES  FOR  TYPE  =  2,  NNELE  =  4 


LN 

SX(I) 

SY(I) 

SXY(I) 

ELEMENT 

LN 

SX(J) 

SY(J) 

SXY(J) 

1 

LN 

SX(K) 

SY(K) 

SXY(K) 

LN 

SX(L) 

SY(L) 

SXY(L) 

ELEMENT 

LN 

SX(I) 

SY(I) 

SXY(I) 

NEL 

LN 

SX(J) 

SY(J) 

SXY(J) 

LN 

SX(K) 

SY  (K) 

SXY(K) 

LN 

SX(L) 

SY  (L) 

SXY(L) 

Figure  3.  Stress  data  file  format 
(Continued) 


8 


STRESSES  FOR  TYPE  =  2.  NNELE 


ELEMENT 

1 


ELEMENT 

NEL 


STRESSES  FOR  TYPE  =  3 


Where  TYPE  = 


ELEMENT 

LN 

SX(M2) 

SY(M2) 

SXY (M2) 

1 

LN 

SX(M3) 

SY(M3) 

SXY (M3) 

LN 

SX(M4) 

SY(M4) 

SXY(M4) 

LN 

SX(C) 

SY(C) 

SXY(C) 

LN 

SX(M1) 

SY(M1) 

SXY(Ml) 

ELEMENT 

LN 

SX(M2) 

SY(M2) 

SXY (M2) 

NEL 

LN 

SX(M3) 

SY(M3) 

SXY (M3) 

LN 

SX(M4) 

SY(M4) 

SXY(M4) 

1  for  stresses 

given 

at  element 

centers 

=  1  for  stresses  given  at  nodal  points  (average  values  at  nodes) 

=  2  for  stresses  given  at  nodal  points  for  each  element 

=  3  for  stresses  given  at  element  centers 

and  at  midsides  for  each  element  (4-node  elements  only) 

TC  =  1  positive  stresses  are  tension 
=  0  positive  stresses  are  compression 
SX  =  global  X  stress 
SY  =  global  Y  stress 
SXY  =  global  XY  stress 
NEL  =  number  of  elements 
NNELE  =  number  of  nodes  per  element  (4  or  8) 

NNP  =  number  of  nodal  points 
LN  *  line  number 
C  =  center 

M1,M2,M3,M4  *  midside  positions  halfway  between  respective  corner  nodes 

Figure  3.  (Concluded) 


allows  the  user  to  select  an  individual  material  property  for  the  calculations. 
All  input  values  are  optional  except  for  the  section  number  I.  Thus,  there 
are  four  specific  forms  of  the  command: 

SECTION  I,  XI,  Yl,  X2,  Y2,  M  (Form  1) 

SECTION  I,  XI,  Yl,  X2,  Y2  (Form  2) 

SECTION  I,  M  (Form  3) 

SECTION  I  (Form  4) 

8.  Form  1  requires  the  user  to  enter  six  values  (section  number,  XI 
coordinate,  Yl  coordinate,  X2  coordinate,  Y2  coordinate,  and  material 
property).  The  calculations  are  based  only  on  the  stresses  from  (XI,  Yl)  to 
(X2,  Y2)  with  elements  having  a  material  property  of  M.  (Stresses  from  ele¬ 
ments  whose  material  property  is  not  "M"  are  omitted  in  calculating  the  re¬ 
sults  . ) 

9.  Form  2  requires  the  user  to  enter  five  values.  The  material  prop¬ 
erty  is  omitted,  and  the  calculations  are  based  on  all  stresses  from  (XI,  Yl) 
to  (X2,  Y2) . 

10.  Form  3  requires  the  user  to  enter  two  values  (section  number  and 
material  property).  The  coordinates  of  (XI,  Yl)  and  (X2,  Y2)  are  set  by  the 
user  with  the  cross  hairs  on  the  graphics  terminal.  Again,  the  calculations 
are  based  on  the  stresses  produced  by  elements  having  a  material  property 
which  corresponds  to  M. 

11.  Form  4  requires  the  user  to  enter  only  the  section  number.  The 
coordinates  of  (XI,  Yl)  and  (X2,  Y2)  are  set  by  the  user  with  the  cross  hairs 
on  the  graphics  terminal.  The  shear,  moment,  and  thrust  calculations  are 
based  on  all  stresses  across  the  selected  section. 

12.  If  XI,  Yl,  X2,  and  Y2  are  not  given  with  the  SECTION  command,  the 
cross  hairs  appear  for  the  user  to  first  pick  (XI,  Yl).  After  the  cross  hairs 
appear,  the  thumb  wheels  on  the  graphics  terminal  can  be  used  to  move  the 
cross  hairs  to  the  location  of  (XI,  Yl).  After  the  user  types  any  character 
(and  possibly  a  carriage  return,  depending  on  how  the  terminal  is  set  up),  the 
cross  hairs  reappear  so  this  process  can  be  repeated  for  (X2,  Y2).  The  grid 
or  a  window  of  the  grid  must  first  be  displayed  before  the  cross  hairs  options 
can  be  used.  The  program  reorders  points  (XI,  Yl)  and  (X2,  Y2)  so  that  the 
first  point  (XI,  Yl)  is  lower  than  the  second  point  (X2,  Y2) .  For  a  horizontal 
line,  the  first  point  is  to  the  left  of  the  second  point.  When  the  section  is 
plotted,  the  first  point  is  always  the  lowest  point. 


10 


PLOT 


13.  This  command  has  two  options: 


PLOT  GRID 
PLOT  SECTION  I 

The  first  command  allows  the  user  to  plot  the  FE  grid.  The  second  plots  the 
Nth  section.  As  many  sections  as  are  defined  can  be  plotted.  If  I  is  not 
specified,  all  the  sections  currently  defined  will  be  plotted. 

WINDOW 

14.  This  command  allows  the  user  to  obtain  a  magnified  view  of  a  rec¬ 
tangular  portion  of  the  grid  by  defining  a  window.  After  typing  the  command, 
the  cross  hairs  appear  for  the  user  to  select  the  lower,  left-hand  corner  of 
the  window.  After  selection,  the  user  types  any  character  and  a  carriage  re¬ 
turn.  The  cross  hairs  then  reappear  so  the  process  can  be  repeated  for  the 
upper  right-hand  corner  of  the  window.  Afterward,  the  window  is  then  plotted. 
OUTPUT 

15.  The  command  is  for  displaying  the  normal  stress,  thrust,  bending 
stress,  and  shear  stress  on  one  plot.  This  command  is  also  given  to  display 
each  stress  distribution  plot  individually  with  the  maximum  and  minimum  stress 
values  displayed. 

lo.  The  form  for  the  combined  plot  is 
OUTPUT  I 

where  I  is  the  section  number.  This  command  allows  the  user  to  plot  the 
stress  distribution  and  resulting  shear,  moment,  and  thrust  for  the  Nth  speci¬ 
fied  section.  Figure  4  shows  a  sample  plot.  The  user  is  given  the  following 
message  if  an  incorrect  section  number  is  given: 

SECTION  NOT  DEFINED 

17.  The  form  for  individual  stress  plots  is 

OUTPUT  NORMAL  STRESS  I 

THRUST 

BENDING  STRESS 
SHEAR  STRESS 
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(XI,  VI)  .  (32.,S.) 
<X2,  V3)  •  (32., It. ) 
NEUTRAL  AXIS  -  (32., 5. ) 
SHEW)  -  -9995. 
NOHEHT  -  -.614SE46 
THRUST  -  •. 


SECTION  NO.  1 


Figure  A.  Section  output  plot 


where  I  is  the  optional  section  number.  If  i  is  not  given,  the  selected 
stress  distribution  will  be  given  for  the  last  section  displayed.  The  follow¬ 
ing  command  was  used  to  generate  Figure  5: 

OUTPUT  I  1 

END 

18.  This  command  terminates  the  program. 


NORflAl  STRESS 


SECTION  NO.  1 

Figure  5.  Normal  stress  plot 
File  Output 

19.  For  each  section  defined,  two  lines  of  data  are  written  to  a  perma¬ 
nent  disc  file.  The  first  line  of  data  contains  the  section  number  and  the 


(XI,  Yl),  (X2,  Y2),  and  (XN,  YN)  coordinates  representing  where  the  neutral 
axis  intersects  the  section.  The  second  line  contains  the  values  of  shear, 
moment,  and  thrust.  See  Appendix  B  for  an  example  display  of  an  output  file. 

Sataple  Run 

20.  Appendix  A  gives  a  step-by-step  example  on  the  use  of  program  CSMT. 
It  shows  uses  of  every  command  but  does  not  demonstrate  every  form  of  each 
command . 


Caution  in  Use  of  the  Program 

21.  The  obvious  use  for  this  program  is  to  calculate  the  shear,  moment, 
and  thrust  along  a  given  section  based  on  output  from  a  FE  analysis.  The  user 
needs  to  remember  that  the  FE  method  is  an  approximate  procedure  that  first 
computes  displacements  and  then  stresses.  These  stresses  can  be  crude  if  an 
adequate  mesh  is  not  used  in  the  analysis.  CSMT  calculates  shear,  moment,  and 
thrust  along  a  section  based  on  interpolation  of  the  computed  discrete  stress 
values.  Because  of  these  assumptions,  equilibrium  may  not  be  satisfied,  par¬ 
ticularly  when  sections  are  taken  in  areas  of  high  stress  gradients  and  es¬ 
pecially  when  a  coarse  grid  is  used.  Results  from  CSMT  at  points  where  loads 
have  been  applied  will  be  of  little  value. 

22.  As  the  grid  in  areas  with  high  stress  gradients  becomes  finer,  re¬ 
sults  from  CSMT  will  more  closely  satisfy  equilibrium.  In  reality,  CSMT  can 
be  used  in  determining  the  validity  of  a  FE  analysis  by  comparing  computed 
values  with  those  required  for  equilibrium. 

Theory  and  Procedure 

23.  The  technique  used  to  compute  shear,  moment,  and  thrust  along  a 
section  is  as  follows: 

a.  Compute  the  principal  X  and  Y  stresses  and  shearing  stress 

(ct  ,  o  ,  T  )  at  the  nodes  and,  when  needed,  at  the  centers 
v  x  ’  y  ’  xy '  * 

of  the  elements. 

b.  Compute  these  stresses  by  interpolation  at  points  along  the 
given  section  where  it  intersects  the  grid. 


_c.  Transform  the  stresses  to  a  new  coordinate  system  with  the  new 
Y  axis  parallel  to  the  section. 

d. -  Compute  thrust  and  shear  using  numerical  integration. 

e .  Compute  the  bending  stress. 

£.  Compute  the  neutral  axis. 

£.  Compute  the  moment. 

Stresses  at  the  nodes  and  centers 


24.  The  first  step  in  the  computation  process  is  to  determine  single 


values  of  normal  X  and  Y  stresses  and  shear  stress  (o  ,  o  ,  and  I  )  for 

v  x  »  y  *  xy' 

each  node  of  the  grid.  For  8-node  elements,  the  stresses  at  the  centers  of 
the  elements  are  also  required. 

25.  TYPE  =  0,  NNELE  =  4.  For  TYPE  =  0,  the  stresses  are  known  at  the 
center  and  must  be  computed  at  the  nodes.  Each  of  the  three  stresses  are 
computed  in  the  same  way;  the  stress  o  at  node  N  (Figure  6)  is  computed  from 


Figure  6. ,  Sample  grid 


Here,  ,  Y^  ,  and  cr  represent  the  coordinates  and  stress,  respectively, 
at  center  i  .  This  is  the  familiar  inverse  distance  squared  interpolation 
formula.  The  summation  (1  to  4  here)  is  in  general  from  1  to  the  number  of 
elements  surrounding  a  node. 

26.  TYPE  -  0,  NNELE  =  8.  The  procedure  is  the  same  as  above,  except 
that  the  stresses  at  the  midside  nodes  must  also  be  computed. 

27.  TYPE  -  1,  NNELE  =  4.  The  stresses  are  input  for  the  nodes  and  are 
not  needed  for  the  elements,  so  no  computation  is  necessary. 

28.  TYPE  =  1 ,  NNELE  =  8 .  The  stresses  are  input  for  the  nodes,  but  the 
stresses  at  the  centers  need  to  be  computed.  The  isoparametric  element  formu¬ 
lation  is  used  and  at  s  =  0  ,  t  =  0  reduces  to  (see  Figure  7) 

°c  =  1  2  °i  ‘  5  2  ai  (2) 

i=5  i=l 

4  7  3 
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8 

Figure  7.  Grid  system  for  8-node  element 

29.  TYPE  =  2,  NNELE  =  4.  This  type  has  as  input  a  stress  for  each  node 

from  each  element.  All  the  different  values  are  simply  averaged  at  a  given 
node  to  determine  a  single  .  Stress  values  at  the  center  are  not  needed. 

30.  TYPE  =  2,  NNELE  =  8.  The  stresses  at  the  nodes  are  computed  by  the 
same  averaging  process  described  for  NNELE  =  4.  The  stress  at  the  center  of 
an  element  is  computed  using  Equation  2  and  the  stress  values  given  for  that 
element. 

31.  TYPE  =  3,  NNELE  =  4.  The  stress  values  at  the  center  and  midside 


positions  are  given  for  each  element.  The  program  treats  the  grid  as  if 
NNELE  =  8.  Thus,  node  numbers  for  the  midside  nodes  are  generated,  and  the  el¬ 
ement  data  are  modified  to  reflect  8  nodes  per  element.  The  stress  values  at 
the  midside  nodes  are  computed  using  a  simple  average  of  the  contribution  from 
each  element.  Equation  1  is  used  to  compute  the  stress  values  at  the  corner 
nodes  with  the  midside  nodes  only  being  used  in  the  summation  (see  Figure  8). 


Figure  8.  Corner  node  computation 
Stresses  along  the  section 

32.  The  next  step  in  the  calculation  procedure  is  to  compute  stress 
values  at  points  along  the  section.  The  points  used  are  those  that  intersect 
the  grid  system.  Figure  9  shows  the  grid  system  used  for  TYPE  =  0,  NNELE  =  4, 
and  Figure  10  shows  the  grid  system  for  TYPE  =1,2,  NNELE  =  4.  For  TYPE  =  3, 
NNELE  =  4,  and  for  all  TYPE  values  for  NNELE  =  8,  the  grid  system  shown  in 
Figure  11  is  used.  The  coordinates  and  stress  values  along  the  section  are 
computed  using  linear  interpolation. 

New  coordinate  system 

33.  Next,  the  coordinate  system  is  rotated  such  that  the  Y'  axis  is 
parallel  to  the  section,  and  the  X'  axis  is  perpendicular  to  the  section  (see 
Figure  12).  The  (X,  Y)  coordinates  are  transformed  by 


X' 


Y' 


cos  A  sin  A 

"x 

-sin  A  cos  A_ 

_Y. 

(3) 
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Figure  11.  Grid  system  for  TYPE  =  3,  NNELE  =  4, 
and  for  all  TYPE  values  for  NNELE  =  8 


Figure  12.  New  coordinate  system 


The  stresses  are  transformed  by 


o'  t' 
x  xy  _ 

t'  o' 

Lxy  y 


cos  A 
-sin  A 


lin  A  CTx  Txy  cos  A 

:os  A  Tv»  sin  A 

j  l  y  j  l 


cos  A  -sin  A 


cos  A 


The  angle  A  is  computed  such  that  its  maximum  absolute  value  is  90  deg.  A 
horizontal  segment  has  A  =  90  deg. 

Spline  fit 


34.  There  now  are  three  sets  of  data  (Y',  o'),  (Y',  o'),  and  (Y',  t  ) 

x  y  xy 

along  the  section.  These  are  first  sorted  and  any  duplicate  points  discarded 


1 sasaesss*  £&$&&&  *&&?&&. 


The  first  and  last  points  are  also  discarded  because  of  end  effects  giving 
incorrect  stress  values  and  are  replaced  using  linear  extrapolation.  Simple 
cubic  splines  are  then  fitted  through  the  stress  points  to  produce  the  three 
diagrams  (for  plotting  purposes  only) . 

Thrust  and  shear 

35.  If  L  is  the  distance  from  the  first  Y'  point  to  the  last  Y' 
point  on  the  section,  then  the  thrust  is  computed  using  numerical  integration 


Y'  +  L 


o'  dY' 
x 


Hence  Yj  is  the  smallest  Y'  value. 

36.  Shear  is  computed  the  same  way  from 


Yi  +  L 


X’  dY' 
xy 


Bending  stress  and  neutral  axis 

37.  The  bending  stress  is  simply 

°b  =  -  r  (7) 

The  neutral  axis  is  where  the  bending  stress  curve  intersects  the  Y'  axis 
and  where  an  equal  area  of  positive  stress  is  separated  from  the  negative  stress 
area.  Call  this  value  Y^  .  X^  is  the  perpendicular  distance  from  the  origin 
of  the  grid  to  the  section.  The  transformation  back  to  the  X-Y  plane  is 

fX„l  f cos  A  -sin  Al  fXil 


sin  A  cos  A 


Moment 


38.  The  moment  is  computed  (using  numerical  integration)  by  integrating 
the  bending  stress  times  moment  arm  along  the  length  L  : 


(9) 


I 


I 


Yj  +  L 

M  =  (Y£  -  Y’)  oB  dY* 

Y' 

X1 

Verification  of  the  Program 

39.  The  cantilever  beam  shown  in  Appendix  A  can  be  used  to  verify  the 
results  from  CSMT.  The  moment  is  linear  in  this  case  and  can  be  computed  from 

M  =  PX 

where 

M  =  moment 

P  =  load  at  the  end  of  the  beam 

X  =  distance  from  the  load  to  the  point  where  the  moment  is  to  be  computed 
The  total  shear  is  constant  at  each  cross  section  (only  vertical  shear  is 
considered)  and  is  equal  to  the  load  at  the  end  of  the  beam.  The  shear 
distribution  is  parabolic.  The  thrust  is  zero  at  all  cross  sections. 

40.  If  sections  are  taken  10,  14,  and  48  ft*  from  the  end  of  the  beam 
as  shown  in  Figure  14,  the  moments  are  calculated  to  be  100,  140,  and 

480  kip-ft.  Figures  15-17  show  the  moments  as  computed  by  CSMT  to  be  96.33, 
135.0,  and  463.2  kip-ft.  These  small  errors  are  a  result  of  the  finite  element 
discretization  of  the  problem  and  interpolation  of  data  by  the  CSMT  program. 

41.  Table  1  shows  a  summary  of  shears  and  moments  with  percent  errors. 

The  calculated  shears,  moments,  and  thrusts  are  for  sections  taken  at  10,  14, 
and  48  ft. 


0 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 
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THRUST 


(ENDING  STRESS 


SHEARING  STRESS 


(HI.  VI)  •  (86  .-I  ) 

(X2.  V2)  .  186. 18  > 

NEUTRAL  AXIS  .  (86  .4  999) 

SHEAR  .  -  1099£*S 
NONENT  •  -  13S9E*6 
THRUST  •  -13  5 

SECTION  NO  a 

CSMT  results  for  section  14  ft  from  beam  end 


THRUST 


SHEARING  STRESS 


(XI.  VI)  •  (Sa  .-1  ) 

<xa.  va>  •  (sa  .is  ) 

NEUTRAL  AXIS  •  (S3  .5  ) 
SHEAR  •  -  1993E«S 
NOHCNT  •  -  4633*6 
THRUST  •  9619 

SECTION  NO  3 


CSMT  results  for  section  48  ft  from  beam  end 


Shear 

Moment 

Distance  from 
Load,  ft 

Calculated 

lb 

Computed 
by  CSMT 
lb 

Percent 

Error 

Calculated 

lb-ft 

Computed 
by  CSMT 
lb-ft 

Percent 

Error 

10 

-10,000 

-9996. 

0.04 

100,000 

96330. 

3.67 

14 

-10,000 

-10000. 

0.00 

140,000 

135000. 

3.57 

48 

-10,000 

-10030. 

0.30 

480,000 

463200. 

3.50 

1.  The  output  file  from  a  typical  FE  analysis  run  (pages  A2-A9)  for  the 
cantilever  beam  shown  in  Figure  A1  contains  both  the  geometry  and  stress  data. 
This  output  will  be  in  a  different  format  for  each  FE  program.  The  user  must 
then  use  the  editor  available  on  the  computer  where  the  output  file  exists  to 
remove  the  unnecessary  data. 


Figure  Al.  Example  cantilever  beam 

2.  Typical  editing  sessions  for  the  geometry  file  on  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  and  Boeing  Computer  Services 
systems  are  shown  on  page  A10.  Editors  on  different  computers  will  have 
different  commands  but  should  have  the  basic  capabilities  illustrated  here. 
All  data  except  for  the  node  and  element  data  have  been  removed  from  the  file 
shown  on  pages  A11-A13.  This  file  (GE0M1)  now  contains  all  node  and  element 
data  but  needs  to  be  rewritten  in  order  to  be  read  by  CSMT.  The  simple  pro¬ 
gram  for  the  WES  computer  shown  on  page  A14  will  generate  the  final  geometry 
data  file  (FGE0M1)  shown  on  pages  A15-A17 . 

3.  Typical  editing  sessions  on  the  WES  and  Boeing  systems  for 
separating  the  stress  data  from  the  output  file  are  listed  on  page  A18.  The 
stress  data  file  (STRF1)  which  was  created  by  editing  the  output  file  (pages 
A2-A9)  is  shown  on  pages  A19-A21.  The  file  STRF1  must  now  be  written  for 
the  program  CSMT  to  read.  The  program  for  the  WES  computer  shown  on  page  A22 
will  generate  the  final  stress  data  file  (FSTRF1)  shown  on  pages  A23  and  A24 . 
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OutpuC  File  fron  a  Typical  FE  Run 


BEAM  PROBLEM  -  LITTLE  GRID 


NUMBER  OF  NODAL  POINTS - 1 

NUMBER  OF  ELEMENTS - 1 

NUMBER  OF  DIFF.  MATERIALS--- 

NUMBER  OF  PRESSURE  CARDS - 

X- ACCELERATION - 

y-acceleration - 

REFERENCE  TEMPERATURE - 

NUMBER  OF  APPROXIMATIONS - 


MAT.  NO.  1  NO.  OF  TEMP.  CARDS  *  1  MASS  DENSITY 

TEMP.  E<C>  NU  E(T)  G/H2 

0,  0.432E  09  0.200E  00  0.432E  09  0.  0. 


STRESS 


DE 

TYPE 

X 

Y 

1 

l  .0 

0. 

10.0 

2 

1.0 

0. 

7.5 

3 

0. 

4.0 

10.0 

4 

3.0 

0. 

5.0 

5 

0. 

4.0 

7.3 

6 

0. 

8.0 

10.0 

7 

1.0 

0. 

2.5 

3 

0. 

4.0 

3.0 

9 

0. 

8.0 

7.3 

10 

0. 

12.0 

10.0 

11 

1.0 

0. 

0. 

12 

0. 

4.0 

2.3 

13 

0. 

8.0 

5.0 

1. 

0. 

12.0 

7.3 

13 

0. 

16.0 

10.0 

16 

0. 

4.0 

0. 

17 

0. 

8.0 

2.3 

18 

0. 

12.0 

3.0 

19 

0. 

16.0 

7.5 

20 

0. 

20.0 

10.0 

21 

0. 

8.0 

0. 

22 

0. 

12.0 

2.5 

23 

0. 

16.0 

3.0 

24 

0. 

20.0 

7.3 

23 

0. 

24.0 

10.0 

26 

0. 

12.0 

0. 

27 

0. 

16.0 

2.5 

28 

0. 

20.0 

3.0 

29 

0. 

24.0 

7.5 

30 

0. 

28.0 

10.0 

31 

0. 

16.0 

0. 

32 

0. 

20.0 

2.5 

33 

0. 

24.0 

3.0 

34 

0. 

28.0 

7.3 

33 

0. 

32.0 

10.0 

36 

0. 

20.0 

0. 

37 

0. 

24.0 

2.5 

38 

0. 

28.0 

3.0 

39 

0. 

32.0 

7.3 

40 

0. 

36.0 

10.0 

41 

0. 

24.0 

0. 

42 

0. 

28.0 

2.3 

43 

0. 

32.0 

3.0 

44 

0. 

36.0 

7.3 

45 

0. 

40.0 

10.0 

46 

0. 

28.0 

0. 

47 

0. 

32.0 

2.5 

48 

0. 

36.0 

3.0 

49 

0. 

40.0 

7.3 

30 

0. 

44.0 

10.0 

51 

0. 

32.0 

0. 

32 

0. 

36.0 

2.3 

33 

0. 

40.0 

3.0 

34 

0. 

44.0 

7.3 

53 

0. 

48.0 

10.0 

36 

0. 

36.0 

0. 

37 

0. 

40.0 

2.3 

38 

0. 

44.0 

3.0 

X-LOAD/DISP 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


Y-LOAD/DISP 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0* 

0. 

0* 

0. 

0. 

0. 

0. 
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59 

0? 

48.0 

7,5 

0. 

0. 

0. 

60 

0. 

52.0 

10.0 

0. 

0. 

0. 

tf 

61 

0. 

40.0 

0. 

0. 

0. 

0. 

r 

62 

0. 

44.0 

2.5 

0. 

0. 

0. 

k 

63 

0. 

48.0 

5.0 

0. 

0. 

0. 

L 

64 

0. 

32.0 

7.5 

0. 

0. 

0. 

65 

0. 

56.0 

10.0 

0. 

0. 

0. 

w 

66 

0. 

44.0 

0. 

0. 

0. 

0. 

2 

67 

0. 

48.0 

2.5 

0. 

0. 

0. 

1 

68 

0. 

52.0 

3.0 

0. 

0. 

0. 

69 

0. 

56.0 

7.5 

0. 

0. 

0. 

Jl 

70 

0. 

60.0 

10.0 

0. 

0. 

0. 

9 

71 

0. 

48.0 

0. 

0. 

0. 

0. 

I 

72 

0. 

52.0 

2.5 

0. 

0. 

0* 

73 

0. 

56.0 

5.0 

0. 

0. 

0. 

74 

0. 

60.0 

7.5 

0. 

0. 

0. 

75 

0? 

64.0 

10.0 

0. 

0. 

0. 

t 

76 

0. 

52.0 

0. 

0. 

0. 

0. 

77 

0. 

56.0 

2.5 

0. 

0. 

0. 

• 

78 

0. 

60.0 

5.0 

0. 

0. 

0. 

79 

0. 

64.0 

7.5 

0. 

0. 

0. 

\ 

80 

0. 

68.0 

10.0 

0. 

0. 

0. 

r 

81 

0. 

56.0 

0. 

0. 

0. 

0. 

82 

0. 

60.0 

2.5 

0. 

0. 

0. 

83 

0. 

64.0 

5.0 

0. 

0. 

0. 

84 

0. 

68.0 

7.5 

0. 

0. 

0. 

85 

0* 

72.0 

10.0 

0. 

0. 

0. 

1 

86 

0. 

60.0 

0. 

0. 

0. 

0. 

87 

0. 

64.0 

2.5 

0. 

0. 

0. 

r 

88 

0. 

68.0 

5.0 

0. 

0. 

0. 

K 

89 

0. 

72.0 

7.5 

0. 

0. 

0. 

90 

0. 

76.0 

10.0 

0. 

0. 

0. 

1 

91 

0. 

64.0 

0. 

0. 

0. 

0. 

92 

0. 

68.0 

2.5 

0. 

0. 

0. 

93 

0. 

72.0 

5.0 

0. 

0. 

0. 

94 

0* 

76.0 

7.5 

0. 

0. 

0. 

95 

0. 

80.0 

10.0 

0. 

0. 

0. 

i 

76 

0. 

68.0 

0. 

0. 

0. 

0. 

97 

0. 

72.0 

2.5 

0. 

0. 

0. 

98 

0. 

76.0 

5.0 

0. 

0. 

0. 

99 

0. 

80.0 

7.5 

0. 

0. 

0. 

100 

0. 

84.0 

10.0 

0. 

0. 

0. 

| 

101 

0. 

72.0 

0. 

0. 

0. 

0. 

102 

0. 

76.0 

2.5 

0. 

0. 

0. 

103 

0. 

80.0 

5.0 

0. 

0. 

0. 

f 

104 

0. 

84.0 

7.3 

0. 

0. 

0. 

103 

0. 

88.0 

10.0 

0. 

0. 

0. 

i 

106 

0. 

76.0 

0. 

0. 

0. 

0. 

107 

0. 

80.0 

2.5 

0. 

0. 

0. 

108 

0. 

84.0 

5.0 

0. 

0. 

0. 

109 

0. 

88.0 

7.3 

0. 

0. 

0. 

f 

110 

0. 

?2 . 0 

10.0 

0. 

0. 

0. 

% 

111 

0. 

80.0 

0. 

0. 

0. 

0. 

112 

0. 

84.0 

2.5 

0. 

0. 

0. 

113 

0. 

88.0 

5.0 

0. 

0. 

0. 

114 

0. 

92.0 

7.5 

0. 

0. 

0. 

115 

0. 

96.0 

10.0 

0. 

0, 

0. 

116 

0. 

84.0 

0. 

0. 

0. 

0. 

) 

117 

0. 

88.0 

2.5 

0. 

0. 

0. 

118 

0. 

?2.0 

5.0 

0. 

0. 

0. 

n? 

0. 

96.0 

7.3 

0. 

0. 

0. 

i 

120 

0. 

100.0 

10.0 

0. 

-0.12SE  04 

0. 

i 

121 

0. 

88.0 

0. 

0. 

0. 

0. 

> 

122 

0. 

92.0 

2.5 

0. 

0. 

0. 

123 

0. 

96.0 

3.0 

0. 

0. 

0. 

124 

0. 

100.0 

7.3 

0. 

-0.250E  04 

0. 

125 

0. 

92.0 

0. 

0. 

0. 

0. 

126 

0. 

96.0 

2.5 

0. 

0. 

0. 

127 

0* 

100.0 

5.0 

0. 

-0.250E  04 

0. 

128 

0. 

96.0 

0. 

0. 

0. 

0. 

12? 

0. 

100.0 

2.5 

0. 

-0.250E  04 

0. 

130 

0. 

100.0 

0. 

0. 

-0.125E  04 

0. 
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EL.  NO. 

I 

J 

K 

L 

HAT. 

W 

1 

11 

16 

12 

7 

1 
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16 

21 

17 
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1 

? 

3 

21 

26 

22 

17 

1 

r 
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26 

31 
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22 
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3 

31 

36 
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27 

t 

fl 

4 

36 

41 
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32 

1 
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46 

42 

37 

1 

l 
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46 

51 

47 

42 

1 

1 
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51 

56 

52 

47 

1 

P 

10 

36 

61 

57 

52 

1 

|i 

l 

A3 

1 

‘  r*  *  -  *  ‘ 

V-  V  V  .• 

NODE 

X-DISP. 

Y-DISP. 

1 

0. 

-0.7832E-04 

2 

0. 

-0 » 2853E-04 

3 

0.4949E-03 

-0 . 2909E-03 

4 

0. 

0. 

5 

0.2447E-03 

-0.2459E-03 

6 

0 . 9A78E-03 

-0.8972E-03 

7 

0. 

-0.2835E-04 

8 

-0.1187E-09 

-0.2331E-03 

9 

0.4794E-03 

-0.8541E-03 

10 

0 • 1 420E-02 

-0.1873E-02 

11 

0. 

-0.7B32E-04 

12 

-0.2447E-03 

-0.2439E-03 

13 

-0.2269E-09 

-0.8395E-03 

14 

0.7055E-03 

-0.1B31E-02 

15 

0  » 1 851 E -02 

-0.3202E-02 

16 

-0.4949E-03 

-0.2909E-03 

17 

-0.4794E-03 

-0.8541E-03 

18 

-0.3362E-09 

-0.1818E-02 

19 

0.9214E-03 

-0 . 3162E-02 

20 

0.2263E-02 

-0.486BE-02 

21 

-0.9678E-03 

-0.8972E-03 

22 

-0.7035E-03 

-0.1831E-02 

23 

-0.4632E-09 

-0.3149E-02 

24 

0.1127E-02 

-0.4831E-02 

23 

0*  2655E-02 

-0.6836E-02 

26 

-0  •  1 420E-02 

-0.1873E-02 

27 

-0.9214E-03 

-0.3162E-02 

28 

-0.5981E-09 

-0. 4818E-02 

29 

0.1323E-02 

-0.6B20E-02 

30 

0.3026E-02 

-0.9149E-02 

31 

-0.18S1E-02 

-0.3202E-02 

32 

-0.1127E-02 

-0.4831E-02 

33 

-0.6927E-09 

-0 • 6808E-02 

34 

0.1S09E-02 

-0.9115E-02 

33 

0.3378E-02 

-0. 1173E-01 

36 

-0.2263E-02 

-0 • 4868E-02 

37 

-0.1323E-02 

-0  *  6820E-02 

38 

-0.8245E-09 

-0.9104E-02 

39 

0.16B3E-02 

-0.1170E-01 

40 

0.3709E-02 

-0. 1439E-01 

41 

-0 . 2655E-02 

-0 *68568-02 

42 

-0.1S09E-02 

-0.91 15E-02 

43 

-0 • 1 131E-08 

-0*1 169E-01 

44 

0. 1830E-02 

-0.1456E-01 

45 

0.4021E-02 

-0.1770E-01 

46 

-0.3026E-02 

-0.9149E-02 

47 

-0.16B5E-02 

-0.U70E-01 

48 

-0.1103E-08 

-0.1433E-01 

49 

0. 2006E-02 

-0.1767E-01 

30 

0.4312E-02 

-0.210SE-01 

51 

-0 . 3378E-02 

-0 . 1 173E-01 

52 

-0.1830E-02 

-0* 1436E-01 

33 

-0.8416E-09 

-0.1766E-01 

34 

0.2132E-02 

-0.2103E-01 

55 

0.43B3E-02 

-0.2463E-01 

36 

-0.3709E-02 

-0* 1459E-01 

57 

-0.2006E-02 

-0.1767E-01 

58 

-0.9294E-09 

-0.2102E-01 

59 

0.2287E-02 

-0.2461E-01 

60 

0.4834E-02 

-0.2842E-01 

61 

-0.4021E-02 

-0.1770E-01 

62 

-0.2152E-02 

-0.2103E-01 

63 

-0.1480E-08 

-0.2460E-01 

64 

0.2413E-02 

-0.2B40E-01 

65 

0.5063E-02 

-0.3240E-01 

66 

-0.4312E-02 

-0.2103E-01 

67 

-0.2287E-02 

-0.2461E-01 

68 

-0.1471E-08 

-0.2839E-01 

69 

0.2328E-02 

-0.323BE-01 

70 

0.3276E-02 

-0*  36S6E-01 

71 

-0.4383E-02 

-0.2463E-01 

72 

-0.2413E-02 

-0.2840E-01 

73 

-0.1208E-08 

-0.3237E-01 

74 

0.2634E-02 

-0.3634E-01 

73 

0.3467E-02 

-0 « 4087E-01 

76 

-0.4834E-02 

-0.2842E-01 

77 

-0.2S28E-02 

-0 ■ 3238E-01 

'..1 


k»' 

W.  *  *. 


I 


A5 


v. 

i 


7B 

-o.  1 1  uc  os 

-0.34531-01 

79 

0.272*1-02 

0.40841-01 

ao 

o.s»m  0? 

0.45341-01 

ii  i 

-0.5063E-02 

-0.32401-01 

02 

-0,24341-0? 

-0.34541-01 

83 

-C.I4IIC-0S 

-0.40851-01 

04 

0.28131-02 

-0.45321-01 

85 

0.578*1-02 

-0.4**31-01 

06 

-0.52761-02 

-0.3*541-01 

87 

-0.272*1-02 

-0.40861-01 

86 

-O.1068E-0B 

-0.45321-01 

89 

0.28*01-02 

-0.4**21-01 

90 

0.5*1*1-02 

-0.54631-01 

91 

-0.54*71-02 

-0.40871-01 

92 

-0.28151-02 

-0.45321-01 

93 

0.7*241-10 

-0.4**11-01 

94 

0.2*551-02 

-0. 54621-01 

95 

0. *0101-02 

-0.5*431-01 

96 

-0.56381-02 

-0.45341-01 

97 

-0.28*01-02 

-0.4**21-01 

98 

0.23741-0* 

-0.54421-01 

99 

0.30111-02 

-0.5*421-01 

too 

0. *1201-02 

-0.64311-01 

101 

-0.578*1-02 

-0.4**31-01 

102 

-0.2*551-02 

-0.54621-01 

103 

0.17501-09 

-0.5*421-01 

104 

0.30561-02 

-0.64311-01 

105 

0. 41911-02 

-0.6*261-01 

106 

-0.591*1-02 

-0.34631-01 

107 

-0.30111-02 

-0.5*421-01 

108 

0.25501-0* 

-0.64301-01 

109 

0.30*11-02 

-0.4*251-01 

no 

0.42411-02 

-0.74251-01 

in 

-0.40301-02 

-0.5943E-01 

112 

-0.30561-02 

-0.6431E-01 

113 

-0.21021-0* 

-0.4*251-01 

114 

0.31161-02 

-0.74251-01 

115 

0.42711-02 

-0.79201-01 

116 

-0.61201-02 

-0.44311-01 

117 

-0.30*11-02 

-0.49251-01 

118 

-0.5393E-0* 

-0.74251-01 

119 

0.31311-02 

-0.7*201-01 

120 

0.42821-02 

-0.04321-01 

121 

-0.41*11-02 

-0.69261-01 

122 

-0.311*1-02 

-0.74231-01 

123 

-0. 74**1-09 

-0.79201-01 

124 

0.313*1-02 

-0.04321-01 

125 

-0.62411-02 

-0. 74251-01 

126 

-0.31311-02 

-0.7*201-01 

127 

0.14861-09 

-O.0432E-O1 

128 

-0.42711-02 

-0.7*201-01 

129 

-0.31361-02 

-0.0432E-01 

130 

-0.62821-02 

-0.0432E-01 

EL  * 

X 

X-STRE0S 

max-stress 

XY-0TREBS 

NO. 

Y 

T-STRESB 

MIN-STREBB 

AN0LE 

1 

2.0 

-0.3**E  03 

0.2171  03 

-0.66*1  03 

1*3 

0.206E  03 

-0.3*91  05 

-0.B92E  02 

2 

6*0 

-0.3021  03 

-0.10*1  02 

-0.7361  03 

1.3 

-0.3311  02 

-O.302E  05 

-O.091E  02 

3 

10*0 

-0.3661  05 

0.20 IE  00 

-0.6531  03 

1*3 

-0.11SE  02 

-0.366E  05 

-0.0*11  02 

4 

14.0 

-0.3S0E  05 

0.1S6E  02 

-0.65*1  03 

1*3 

0.3171  01 

-0.350E  05 

-0 . 0*11  02 

5 

18*0 

-0.3331  05 

0.1261  02 

-0.63*1  03 

1*3 

-0.4341  00 

-0.3341  05 

-0.0*01  02 

6 

22.0 

-0.317E  03 

0.137E  02 

-0.6591  03 

1 » 3 

0.6051-01 

-0.317E  05 

-0.0*01  02 

7 

26*0 

-0.3011  05 

0.144E  02 

-0.65*1  03 

1*3 

0.4091-02 

-0.3011  05 

-0.S89E  02 

8 

30*0 

-0.2051  03 

0.1521  02 

-0.65*1  03 

1.3 

-0.3*61-02 

-0.2051  OS 

-0.888E  02 

9 

34*0 

-0.2601  03 

0.162E  02 

-0.639E  03 

1*3 

-0.134E-01 

-0.26*1  05 

-0.8071  02 

10 

38.0 

-0.232E  05 

0.172E  02 

-0.659E  03 

1.3 

-0.3131-01 

-0.2321  05 

-O.807E  02 

A6 


0.23AE  OS 
0.105E-01 
0.220E  05 
0.110E  00 
0.203E  OS 
0  *  453E-01 
0.187E  OS 
0.789E-01 
0.171E  05 
0.73SE-01 
0.1SSE  OS 
0 • 158E  00 
0.138E  OS 
0.131E  00 
0.122E  03 
0.170E  00 
0.106E  OS 
0.134E-01 
0.896E  04 
0.439E-01 
0.733E  04 
0.192E  00 
0.370E  04 
0.988E-01 
0.407E  04 
0.383E  00 
0.244E  04 
0.249E  02 
0.814E  03 
0.172E  03 


X-STRESS 
Y-STRESS 
0.130E  OS 
0.976E  03 
0.127E  OS 
0.170E  03 
0.122E  OS 
0.103E  02 
0.U7E  OS 
0.417E  00 
0.U1E  OS 
0.216E  00 
0.106E  OS 
0.443E-01 
0.100E  OS 
0.321E-01 
0.949E  04 
0.130E-01 
0.89SE  04 
0.212E-01 
0 • 84  IE  04 
0.4S8E-01 
0.786E  04 
0.186E-01 
0.732E  04 
0.121E-03 
0.678E  04 
0.1 33E  00 
0.624E  04 
0.181E  00 
0.S70E  04 
0.107E  00 
0.S1SE  04 
0.763E-01 
0 • 46 IE  04 
0.136E  00 
0.407E  04 
0.364E  00 
0.3S3E  04 
0.803E-01 
0.299E  04 
0.464E-01 
0.244E  04 
0.722E-01 
0.190E  04 
0.431E  00 
0.136E  04 
0.104E  01 
0.803E  03 
0.195E  02 
0.245E  03 


0.184E  02 
-0.236E  05 
0.199E  02 
-0.220E  05 
0.214E  02 
-0.204E  05 
0.231E  02 
-0.187E  OS 
0.2S3E  02 
-0.171E  OS 
0.282E  02 
-0.153E  OS 
0.31SE  02 
-0.139E  05 
0.3S4E  02 
-0.122E  05 
0.410E  02 
-0.106E  05 
0.483E  02 
-0.900E  04 
0.588E  02 
-0.739E  04 
0.7SSE  02 
-0.578E  04 
0.105E  03 
-0.418E  04 
0.135E  03 
-0.260E  04 
0.S17E  03 
-0.116E  04 


MAX-STRESS 
HEN-STRESS 
0.110E  04 
-0.131E  OS 
-0.448E  02 
-0.1 28E  OS 
0.1S6E  03 
-0.123E  OS 
0.131E  03 
-0.1 18E  OS 
0.1S9E  03 
-0.113E  OS 
0.166E  03 
-0.107E  05 
0.17SE  03 
-0.102E  OS 
0.18SE  03 
-0.968E  04 
0.19SE  03 
-0.914E  04 
0.207E  03 
-0.861E  04 
0.221E  03 
-0.809E  04 
0.236E  03 
-0.756E  04 
0.254E  03 
-0.703E  04 
0.274E  03 
-0.651E  04 
0.298E  03 
-0.599E  04 
0.326E  03 
-0.548E  04 
0.360E  03 
-0.497E  04 
0.400E  03 
-0.447E  04 
0.450E  03 
-0.398E  04 
0.S12E  03 
-0.3S0E  04 
0.S90E  03 
-0.303E  04 
0.691E  03 
-0.2S9E  04 
0.820E  03 
-0.218E  04 
0.100E  04 
-0.182E  04 
0.128E  04 


-0.659E  03 
-0.886E  02 
-0.659E  03 
-0.B84E  02 
-0.6S9E  03 
-0.883E  02 
-0.6S9E  03 
-0.881E  02 
-0.6S9E  03 
-0.879E  02 
-0.660E  03 
-0.877E  02 
-0.659E  03 
-0.874E  02 
-0.660E  03 
-0.871E  02 
-0.660E  03 
-0.866E  02 
-0.660E  03 
-0.860E  02 
-0.660E  03 
-0.850E  02 
-0.660E  03 
-0.836E  02 
-0.661E  03 
-0.811E  02 
-0.643E  03 
-0.761E  02 
-0.677E  03 
-0.631E  02 


XY-STRESS 
ANGLE 
-0.133E  04 
-0.848E  02 
-0.126E  04 
-0.843E  02 
-0.134E  04 
-0.839E  02 
-0.134E  04 
-0.837E  02 
-0.134E  04 
-0.834E  02 
-0.134E  04 
-0.831E  02 
-0.134E  04 
-0.827E  02 
-0.1 34E  04 
-0.823E  02 
-0.134E  04 
-0.81BE  02 
-0.134E  04 
-0.B13E  02 
-0.134E  04 
-0.808E  02 
-0.134E  04 
-0.801E  02 
-0.134E  04 
-0.794E  02 
-0.134E  04 
-0.783E  02 
-0.134E  04 
-0.776E  02 
-0.134E  04 
-0.764E  02 
-0.134E  04 
-0.751E  02 
-0.1 34E  04 
-0.73SE  02 
-0.134E  04 
-0.713E  02 
-0.134E  04 
-0.692E  02 
-0.134E  04 
-0.663E  02 
-0.134E  04 
-0.62BE  02 
-0.134E  04 
-0.S86E  02 
-0.134E  04 
-0.S31E  02 
-0.132E  04 


EL. 

X 

NO. 

Y 

51 

2.0 

4.3 

52 

4.0 

4.3 

53 

10.0 

4.3 

54 

14.0 

4.3 

55 

18.0 

4.3 

54 

22.0 

4.3 

5? 

24.0 

4.3 

58 

30.0 

4.3 

5? 

34.0 

4.3 

60 

38.0 

4.3 

61 

42.0 

4.3 

62 

44.0 

4.3 

63 

50.0 

4.3 

64 

54.0 

4.3 

65 

58.0 

4.3 

66 

42.0 

4.3 

67 

44.0 

4.3 

68 

70.0 

4.3 

49 

74.0 

4.3 

70 

78.0 

4.3 

71 

82.0 

4.3 

72 

84.0 

4.3 

73 

90.0 

4.3 

74 

94.0 

4.3 

75 

98.0 

4.3 

EL* 

X 

NO. 

Y 

76 

2.0 

8.8 

77 

4.0 

8.8 

78 

10.0 

8.8 

7? 

14.0 

8.8 

80 

18.0 

8.8 

81 

22.0 

8.8 

82 

24.0 

8.8 

83 

30.0 

8.8 

84 

34.0 

8.8 

83 

38.0 

8.8 

86 

42.0 

8.8 

87 

44.0 

8.8 

X-STRESS 
Y-STRESS 
0.130E  05 
-0.974E  03 
0.127E  05 
0.170E  03 
0.122E  05 
-0.103E  02 
0.117F.  05 
0.421E  CO 
0.111E  05 
-0.217E  00 
0* 106E  05 
0.435E-01 
O.tOOE  05 
0.173E-01 
0 . 949E  04 
0.502E-01 
0.89SE  04 
-0.421E-01 
0.841E  04 
-0.991E-01 
0.7B4E  04 
0.200E-01 
0.732E  04 
0.994E-01 
0.478E  04 
-0.585E-01 
0.424E  04 
-0.110E  00 
0.570E  04 
-0.289E-02 
0.515E  04 
0.404E-02 
0.441E  04 
-0.202E  00 
0.407E  04 
-0.240E  00 
0.3S3E  04 
-0.582E-01 
0.299E  04 
-0.593E-01 
0.244E  04 
0.428E-01 
0.190E  04 
-0.411E  00 
0.134E  04 
0.980E  00 
0.B03E  03 
0.194E  02 
0.245E  03 
-0.I40E  03 


X-STRESS 
Y-STRESS 
0.399E  05 
-0.20AE  03 
0.382E  05 
0.331E  02 
0.344E  03 
0.113E  02 
0.330E  05 
-0.317E  01 
0.333E  05 
0.434E  00 
0.317E  05 
-0.472E-01 
0.301E  05 
0.223E-01 
0.285E  05 
0.129E-01 
0.248E  05 
0.540E-01 
0.252E  03 
-0.170E-01 
0.23&E  OS 
-0.344E-01 
0.220E  03 
0.149E-01 


MAX-STRESS 
MIN-STRESS 
0.131E  05 
-O.tlOE  04 
0.128E  OS 
0.448E  02 
0.123E  05 
-0.154E  03 
0.118E  05 
-0.131E  03 
0.113E  05 
-0.139E  03 
0.107E  05 
-0.1 A4E  03 
0.102E  05 
-0.173E  03 
0 . 968E  04 
-0.185E  03 
0.914E  04 
-0.195E  03 
0.841E  04 
-0.208E  03 
0.809E  04 
-0.221E  03 
0 . 756E  04 
-0.234E  03 
0.703E  04 
-0.234E  03 
0.451E  04 
-0.275E  03 
0.599E  04 
-0.298E  03 
0.S48E  04 
-0.324E  03 
0.497E  04 
-0.340E  03 
0.447E  04 
-0.400E  03 
0.398E  04 
-0.450E  03 
0.330E  04 
-0.S12E  03 
0.303E  04 
-0.390E  03 
0.239E  04 
-0.491E  03 
0.21BE  04 
-0.820E  03 
0.1B2E  04 
-0.100E  04 
0.139E  04 
-0.128E  04 


MAX-STRESS 
MIN-STRESS 
0.399E  05 
-0.217E  03 
0.382E  05 
0.189E  02 
0.344E  03 
-0.199E  00 
0.350E  05 
-0.154E  02 
0.334E  03 
-0.12AE  02 
0.317E  05 
-0.137E  02 
0.301E  OS 
-0.144E  02 
0.28SE  05 
-0.132E  02 
0.249E  05 
-0.141E  02 
0.232E  05 
-0.172E  02 
0.234E  05 
-0.184E  02 
0.220E  05 
-0.197E  02 


XY-STRESS 
ANGLE 
-0.133E  04 
-0.S38E  01 
-0.124E  04 
-0.S49E  01 
-0.134E  04 
-0.421E  01 
-0.134E  04 
-0.447E  01 
-0.134E  04 
-0.477E  01 
-0.134E  04 
-0.711E  01 
-0.134E  04 
-0.747E  01 
-0.134E  04 
-0.78BE  01 
-0.134E  04 
-0.833E  01 
-0.134E  04 
-0.884E  01 
-0.134E  04 
-0.940E  01 
-0.134E  04 
-0.100E  02 
-0.134E  04 
-0.108E  02 
-0.134E  04 
-0.114E  02 
-0.134E  04 
-0.124E  02 
-0.134E  04 
-0.1 37E  02 
-0.1 34E  04 
-0.131E  02 
-0.1 34E  04 
-0.147E  02 
-0.134E  04 
-0.184E  02 
-0.134E  04 
-0.210E  02 
-0.134E  04 
-0.238E  02 
-0.134E  04 
-0.274E  02 
-0.134E  04 
-0.316E  02 
-0.134E  04 
-0.370E  02 
-0.1 32E  04 
-0.409E  02 


XY-8TRES8 
ANGLE 
-0.449E  03 
-0.9S7E  00 
-0.734E  03 
-0.111E  01 
-0.433E  03 
-0.102E  01 
-0.4S9E  03 
-0.108E  01 
-0.4S9E  03 
-0.113E  01 
-0.459E  03 
-0.H9E  01 
-0.459E  03 
-0.124E  01 
-0.459E  03 
-0.133E  01 
-0.459E  03 
-0.141E  01 
-0.459E  03 
-0.150E  01 
-0.459E  03 
-0.140E  01 
-0.459E  03 
-0.172E  01 
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Listing  of  Program  on  WES  Computer  for  Generatin 
the  Geometry  Data  File 


1000C 

1010C  PROGRAM  FOR  CREATING  GEOMETRY  DATA  FILE. 

1020C 

1030  DIMENSION  IE(6> 

1040  CHARACTER  NAME*6,  FNAME48 

1050  LIN  -  10000 

1060C 

1070C  ATTACH  DATA  FILE  CONTAINING  GEOMETRY  DATA. 

1080C 

1090  PRINT.  'ENTER  INPUT  FILENAME* 

1100  READ  100.  NAME 

1110  100  FORMAT  (A6) 

1120  ENCODE  (FNAME.110)  '/*.  NAME.  *  1  * 

1130  110  FORMAT  (Al.  A6,  Al) 

1140  CALL  ATTACH<  01 »  FNAME .  3.  0.  1ST) 

1150C 

1160C  ATTACH  FILE  FOR  GEOMETRY  TO  BE  WRITTEN 

1170C  INTO  IN  PROPER  FORMAT. 

1180C 

1190  PRINT.  "ENTER  OUTPUT  FILENAME* 

1200  READ  100.  NAME 

1210  ENCODE  < FNAME .110)  */*.  NAME.  *»* 

1220  CALL  ATTACH(02.  FNAME.  3.  0.  1ST) 

1230C 

1240  PRINT.  ‘ENTER  NUHBER  OF  NODE  POINTS' 

1250  READ  120.  NNP 

1260  120  FORMAT  (V) 

1270C 

1280  PRINT.  ‘ENTER  NUMBER  OF  ELEMENTS* 

1290  READ  120.  NEL 

1300C 

1310  PRINT.  'ENTER  NUMBER  OF  NODES  PER  ELEMENT* 

1320  READ  120.  NNELE 

1330C 

1340  PRINT.  ‘ENTER  NUMBER  OF  MATERIAL  TYPES* 

1350  READ  120.  NMAT 

1360C 

1370  WRITE  (2.130)  LIN.  NNP.  NEL,  NNELE,  NMAT 

1380  130  FORMAT  (715) 

1390  LIN  -  LIN  +  10 

1400C 

14 IOC  NODE  DATA 

1420C 

1430  DO  150  I  »  1,  NNP 

1440  READ  (1,120)  NP,  BC,  X,  Y 

1450  WRITE  (2.140)  LIN.  NP,  X,  Y 

1460  140  FORMAT  (215,  2F10.1) 

1470  150  LIN  *  LIN  +  10 

1480C 

1490C  ELEMENT  Her A 

1500C 

1510  DO  160  I  *  1,  NEL 

1520  READ  (1,120)  NE,  IE 

1530  WRITE  (2,130)  LIN,  NE,  IE 

1540  160  LIN  -  LIN  4  10 

1550C 
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-  *J8  * 
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1 

ft  10790 
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64.0 

7.5 

'i  loaoo 

80 

68.0 

10.0 

„•  10910 

81 

56.0 

0. 

ft  10820 

82 

60.0 

2.5 

hk  10830 

83 

64.0 

3.0 

1  10840 

84 

68.0 

7.3 

■  10850 

85 

72.0 

10.0 

■;  10860 

86 

60.0 

0. 

v  10870 

87 

64.0 

2.3 

r  10880 

88 

68.0 

5.0 

£  10890 

89 

72.0 

7.5 

.•  10900 

90 

76.0 

10.0 

•  •  10910 

91 

64.0 

0. 

£  10920 

92 

68.0 

2.3 

|  10930 

93 

72.0 

5.0 

|  10940 

94 

76.0 

7.5 

10930 

95 

80.0 

10.0 

S*.  10960 

96 

68.0 

0. 

£*  10970 

97 

72.0 

2.5 

10980 

98 

76.0 

5.0 

10990 

99 

80.0 

7.5 

;•  nooo 

100 

84.0 

10.0 

C*  11010 

ll  11020 

101 

72.0 

0. 

102 

76.0 

2.3 

|  11030 

103 

80.0 

5.0 

^  11040 

104 

84.0 

7.3 

•J  11030 

103 

88.0 

10.0 

S  11060 

106 

76.0 

0. 

11070 

107 

80.0 

2.5 

11080 

108 

84.0 

3.0 

\  11090 

109 

88.0 

7.5 

N  11100 

110 

92.0 

10.0 

U  11110 

111 

80.0 

0. 

a  11120 

112 

84.0 

2.3 

R  11130 

113 

88.0 

5.0 

S  11140 

114 

92.0 

7.3 

K  11130 

113 

96.0 

10.0 

n  moo 

116 

84.0 

0. 

r  11170 

117 

88.0 

2.5 

B  11180 

118 

92.0 

5.0 

P  11190 

119 

96.0 

7.3 

a  11200 

120 

100.0 

10.0 

1  11210 
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88.0 

0. 

1  11220 
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92.0 

2.3 

H  11230 
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3.0 

11240 
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100.0 

7.3 

r*  11230 
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92.0 
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«  11260 
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3  ”270 
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*  11280 

128 

96.0 

0. 

'i  11290 
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27 
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9 
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32 

47 
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I  11410 
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1 
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1 
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1  11370 

27 

12  17 

13 

8 

1 
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28 
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18 

13 

1 
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# 

29 
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23 

18 
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37 
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58 
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38 

67 

72 
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63 
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39 

72 

77 

73 

68 
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40 

77 

82 

78 

73 
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41 

82 

87 

83 

78 

11720 

42 

87 

92 
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83 
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43 

92 

97 

93 

88 
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44 

97 

102 

98 

93 

11750 

45 
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107 

103 

98 

11760 

46 

107 

112 

108 

103 

11770 

47 
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117 

113 

108 

11780 

48 

117 

122 
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113 

11790 

49 

122 

126 
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118 

11800 

50 

126 

129 

127 
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51 

4 

8 

5 

2 
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52 

8 

13 
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5 
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53 
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9 

11840 

54 

18 

23 

19 

14 

11850 

55 

23 

28 

24 

19 

11860 

56 

28 

33 

29 

24 

11870 

57 

33 

38 

34 

29 

11880 

58 

38 

43 

39 

34 

11890 

59 

43 

48 

44 

39 

11900 

60 

48 

S3 

49 

44 

11910 

61 

S3 

58 

54 

49 

11920 

62 

58 

63 

59 

54 

11930 

63 

63 

68 

64 

59 

11940 

64 

68 

73 

69 

64 

11950 

65 

73 

78 

74 

69 

11960 

66 

78 

83 

79 

74 

11970 

67 

83 

88 

84 

79 

11980 

68 

88 

93 

89 

84 

11990 

69 

93 

98 

94 

89 

12000 

70 

98 

103 

99 

94 

12010 

71 

103 

108 

104 

99 

12020 

72 

108 

113 

109 

104 

12030 

73 

113 

118 

114 

109 

12040 

74 

118 

123 

119 

114 

12050 

75 

123 

127 

124 

119 

12060 

76 

2 

5 

3 

1 

12070 

77 

5 

9 

6 

3 

12080 

78 

9 

14 

10 

6 

12090 

79 

14 

19 

15 

10 

12100 

80 

19 

24 

20 

15 

12110 

81 

24 

29 

25 

20 
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82 

29 

34 

30 

25 
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83 

34 

39 

35 

30 
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84 

39 

44 

40 

35 
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85 

44 

49 

45 

40 
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86 

49 

54 

50 

45 
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87 

54 

59 

55 

SO 

12180 

88 

59 

64 

60 

55 
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89 

64 

69 

65 

60 
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90 

69 

74 

70 

65 

12210 

91 

74 

79 

75 

70 

12220 

92 

79 

84 

80 

75 
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93 

84 

89 

85 

80 
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94 

89 

94 

90 

85 

12250 

95 

94 

99 

95 

90 

12260 

96 

99 
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95 
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97 
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109 
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too 

12280 

98 
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114 
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12290 

99 

114 

119 

115 

110 

12300 

100 

119 

124 

120 

115 
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Typical  Editing  Sessions  for  Separating  the  Stress  Data 
from  the  Output  File 


WES  System 

♦EDIT 

-DI26 

-P 

Boeing  System 
OCMEDIT 

E>DE26 

E>P 

NODE 

TYPE  X  Y 

X-LOAD/DISP 

Y-LOAD/DISP 

-DI13S 

-P 

E>DE135 

E>P 

EL.  NO. 

I  J  K 

L  HAT. 

-D 1 107 
-P 

E>DE107 

E>P 

NODE 

X-DISP •  Y-DISP. 

-DI135 

-P 

E>DE135 

E>P 

EL. 

X  X-STRESS 

HAX-STRESS 

XY-STRESr; 

-DI2 

-P 

E>DE2 

E>P 

1 

2.0  -0.400E  05 

0.218E  03 

-0.670E  03 

-FI50 

-DI3 

-P 

E>N50 

E>DE3 

E>P 

NO. 

Y  Y-STRESS 

HIN-STRESS 

ANGLE 

-Dll 
-F 150 
-DM 
-P 

E>DE1 

E>N50 

E>DE4 

E>P 

SI  2.0  0.130E  05 

-FI  50 

-DM 

-P 

76  2.0  0.400E  05 

-F 150 
-DM 

•nd  of  file  -  reauest  executed 

-SAVE  STRF1 
DATA  SAVED-STRF1 
end  of  file 


0.132E  05  -0.133E  04 

E>N50 

E>DE4 

E>P 

0.400E  05  -0.670E  0;> 

E>N50 

E>D4 

2  ti*?& 

E>File  10 

File  STRF2  edited 

and  replaced 


TEHP. 


-DONE 

* 
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Lilting  of  Str«»»  D«t«  Fll«  STK71 


♦LIST  STRF1 


1 

2.0 

-0.399E  03 

0.217E  03 

-0.669E  03 

1.3 

0.206E  03 

-0.399E  05 

-0.B92E  02 

2 

6.0 

-0.382E  03 

-0.189E  02 

-0.736E  03 

1.3 

-0.331E  02 

-0.382E  05 

-0.891E  02 

3 

10.0 

-0.366E  05 

0.20 IE  00 

-0.6S3E  03 

1.3 

-0.115E  02 

-0.366E  03 

-0.891E  02 

4 

14.0 

-0.350E  03 

0.1S6E  02 

-0.6S9E  03 

1.3 

0.317E  01 

-0.350E  05 

-0.891E  02 

S 

18.0 

-0.333E  05 

0.126E  02 

-0.639E  03 

1.3 

-0.434E  00 

-0.334E  05 

-0.B90E  02 

6 

22.0 

-0.317E  05 

0.137E  02 

-0.6S9E  03 

1.3 

0.605E-01 

-0.317E  05 

-0.890E  02 

7 

26.0 

-0.301E  05 

0.144E  02 

-0.6S9E  03 

1.3 

0 • 409E-02 

-0.301E  05 

-0.889E  02 

e 

30.0 

-0.285E  05 

0.132E  02 

-0.639E  03 

1.3 

-0 . 396E-02 

-0.28SE  03 

-0.88BE  02 

9 

34.0 

-0.268E  05 

0.162E  02 

-0.659E  03 

1.3 

-0.136E-01 

-0.269E  03 

-0.B87E  02 

10 

38.0 

-0.252E  05 

0.172E  02 

-0.6S9E  03 

1.3 

-0.315E-01 

-0.252E  03 

-0.BB7E  02 

u 

42.0 

-0.236E  05 

0.184E  02 

-0.6S9E  03 

1.3 

0. 105E-01 

-0.236E  OS 

-0.886E  02 

12 

46.0 

-0.220E  05 

0.199E  02 

-0.6S9E  03 

1.3 

0.110E  00 

-0.220E  OS 

-0.BB4E  02 

13 

30.0 

-0.203E  05 

0.214E  02 

-0.639E  03 

1.3 

0.433E-01 

-0.204E  03 

-0.883E  02 

14 

34.0 

-0.1B7E  05 

0.231E  02 

-0.639E  03 

1.3 

-0.789E-01 

-0.187E  05 

-0.881E  02 

15 

58.0 

-0.171E  05 

0.2S3E  02 

-0.6S9E  03 

1.3 

-0.753E-01 

-0.171E  OS 

-0.B79E  02 

16 

62.0 

-0.155E  05 

0.282E  02 

-0.660E  03 

1.3 

0.158E  00 

-0.1S5E  03 

-0.B77E  02 

17 

66.0 

-0.138E  03 

0.31SE  02 

-0.659E  03 

1.3 

0.131E  00 

-0.1 39E  03 

-0.874E  02 

18 

70.0 

-0.122E  05 

0.354E  02 

-0.660E  03 

1.3 

-0.1 70E  00 

-0.1 22E  03 

-0.871E  02 

19 

74.0 

-0.106E  03 

0.410E  02 

-0.660E  03 

1.3 

0.134E-01 

-0.1 06E  03 

-0.866E  02 

20 

78.0 

-0.896E  04 

0.483E  02 

-0.660E  03 

1.3 

-0.439E-01 

-0.900E  04 

-0.860E  02 

21 

82.0 

-0.733E  04 

0.388E  02 

-0.660E  03 

1.3 

-0.192E  00 

-0.739E  04 

-0.850E  02 

22 

86.0 

-0.570E  04 

0.733E  02 

-0.660E  03 

1.3 

0.988E-01 

-0.578E  04 

-0.B36E  02 

23 

90.0 

-0.407E  04 

0.103E  03 

-0.661E  03 

1.3 

0.383E  00 

-0.418E  04 

-0.811E  02 

24 

94.0 

-0.244E  04 

0.133E  03 

-0.643E  03 

1.3 

-0.249E  02 

-0.260E  04 

-0.761E  02 

23 

98.0 

-0.814E  03 

0.317E  03 

-0.677E  03 

1.3 

0.172E  03 

-0.1 16E  04 

-0.631E  02 

26 

2.0 

-0.1 30E  05 

0.1 10E  04 

-0.133E  04 

3.8 

0.976E  03 

-0.131E  05 

-0.B48E  02 

27 

6.0 

-0.127E  05 

-0.44SE  02 

-0.126E  04 

3.8 

-0.1 70E  03 

-0.1 28E  03 

-0.843E  02 

28 

10.0 

-0.122E  05 

0.156E  03 

-0.134E  04 

3.8 

0.103E  02 

-0.123E  03 

-0.839E  02 

29 

14.0 

-0.117E  03 

0.151E  03 

-0.134E  04 

3.8 

-0.417E  00 

-0.118E  03 

-0.837E  02 

30 

18.0 

-0. 11 IE  03 

0.159E  03 

-0.1 34E  04 

3.8 

0.216E  00 

-0.U3E  05 

-0.834E  02 

31 

22.0 

-0.106E  05 

0.166E  03 

-0.134E  04 

3.8 

-0.443E-01 

-0.107E  OS 

-0.B31E  02 

32 

26.0 

-0.100E  05 

0.175E  03 

-0.134E  04 

3.8 

0.321E-01 

-0.102E  03 

-0.827E  02 

33 

30.0 

-0.949E  04 

0.185E  03 

-0.134E  04 

3.8 

0.130E-01 

-0.96BE  04 

-0.823E  02 

34 

34.0 

-0.893E  04 

0.19SE  03 

-0.134E  04 

3.8 

-0.212E-01 

-0.914E  04 

-0.B18E  02 

33 

38.0 

-0.841E  04 

0.207E  03 

-0.134E  04 

3.8 

-0.438E-01 

-0.B61E  04 

-0.813E  02 

36 

42.0 

-0.786E  04 

0.221E  03 

-0.134E  04 

3.8 

-0 . 186E-01 

-0.809E  04 

-0.808E  02 

37 

46.0 

-0.732E  04 

0.236E  03 

-0.134E  04 

3.8 

-0.121E-03 

-0.736E  04 

-0.801E  02 

38 

30.0 

-0.678E  04 

0.254E  03 

-0.134E  04 

3.8 

-0.1 33E  00 

-0.703E  04 

-0.794E  02 

3? 

34.0 

-0.624E  04 

0.274E  03 

-0.134E  04 

3.8 

-0.181E  00 

-0.631E  04 

-0.7BSE  02 

40 

38.0 

-0.S70E  04 

0.298E  03 

-0.134E  04 

3.8 

-0.107E  00 

-0.399E  04 

-0.776E  02 

•0.S15E  04 
0.743E-01 
0.461E  04 
■0.134E  00 
0.407E  04 
0.364E  00 
0.353E  04 
0.803E-01 
0.299E  04 
0.464E-01 
0.244E  04 
0.722E-01 
0.190E  04 
0.431E  00 
0.136E  04 
0.104E  01 
0.803E  03 
0.195E  02 
0.243E  03 
0.140E  03 
0.130E  OS 
0.976E  03 
0.127E  OS 
0.170E  03 
0.122E  OS 
0.103E  02 
0.U7E  OS 
0.421E  00 
0.111E  OS 
0.217E  00 
0.104E  OS 
0.435E-01 
0.100E  OS 
0.173E-01 
0.949E  04 
0.S02E-01 
0.89SE  04 
0.421E-01 
0.841E  04 
0.991E-01 
0.786E  04 
0.200E-01 
0.732E  04 
0.994E-01 
0.478E  04 
0.S8SE-01 
0.624E  04 
0.110E  00 
0.S70E  04 
0.289E-02 
0.S13E  04 
0.404E-02 
0 . 4A1E  04 
0.202E  00 
0.407E  04 
0.240E  00 
0.353E  04 
0.S82E-01 
0.299E  04 
0.S93E-01 
0.244E  04 
0.428E-01 
0.190E  04 
0.411E  00 
0.13AE  04 
0.980E  00 
0.803E  03 
0.194E  02 
0.24SE  03 
0.140E  03 
0.399E  05 
0.206E  03 
0.382E  OS 
0.331E  02 
0.366E  OS 
0.11SE  02 
0.3S0E  05 
0.317E  01 
0.333E  05 
0.434E  00 
0.317E  OS 
0.472E-01 


0.326E  03 
-0.S48E  04 
0.360E  03 
-0.497E  04 
0.400E  03 
-0.447E  04 
0.450E  03 
-0.398E  04 
0.S12E  03 
-0.350E  04 
0.390E  03 
-0.303E  04 
0 . 691E  03 
-0.259E  04 
0.820E  03 
-0.218E  04 
0.100E  04 
-0.182E  04 
0.128E  04 
-0.139E  04 
0.131E  05 
-0.110E  04 
0.128E  OS 
0.448E  02 
0.123E  OS 
-0.1S6E  03 
0.118E  05 
-0.1S1E  03 
0.113E  03 
-0.159E  03 
0.107E  OS 
-0.16AE  03 
0.102E  OS 
-0.175E  03 
0.948E  04 
-0.18SE  03 
0.914E  04 
-0.19SE  03 
0.841E  04 
-0.208E  03 
0.809E  04 
-0.221E  03 
0.73AE  04 
-0.23AE  03 
0.703E  04 
-0.2S4E  03 
0.A31E  04 
-0.27SE  03 
0.S99E  04 
-0.298E  03 
0.S48E  04 
-0.326E  03 
0.497E  04 
-0.340E  03 
0.447E  04 
-0.400E  03 
0.398E  04 
-0.4SOE  03 
0.3S0E  04 
-0.S12E  03 
0.303E  04 
-0.390E  03 
0.259E  04 
-0.691E  03 
0.218E  04 
-0.820E  03 
0.182E  04 
-0.100E  04 
0.139E  04 
-0.1 28E  04 
0.399E  OS 
-0.217E  03 
0.382E  OS 
0.189E  02 
0.366E  05 
-0.199E  00 
0.3S0E  05 
-0.1 56E  02 
0.334E  03 
-0.126E  02 
0.317E  OS 
-0.137E  02 


-0.134E  04 
-0.764E  02 
-0.134E  04 
-0.7S1E  02 
-0.134E  04 
-0.735E  02 
-0.134E  04 
-0.715E  02 
-0.1 34E  04 
-0.692E  02 
-0.134E  04 
-0 • 663E  02 
-0.134E  04 
-0.628E  02 
-0.1 34E  04 
-0.S86E  02 
-0.136E  04 
-0.S31E  02 
-0.132E  04 
-0.492E  02 
-0.133E  04 
-0.S38E  01 
-0. 126E  04 
-0.S69E  01 
-0.1 34E  04 
-0.621E  01 
-0.134E  04 
-0.A47E  01 
-0.134E  04 
-0.677E  01 
-0.134E  04 
-0.71  IE  01 
-0.1 34E  04 
-0.747E  01 
-0.134E  04 
-0.788E  01 
-0.134E  04 
-0.833E  01 
-0.134E  04 
-0.884E  01 
-0.134E  04 
-0.940E  01 
-0.I34E  04 
-0.100E  02 
-0.134E  04 
-0.108E  02 
-0.134E  04 
-0.116E  02 
-0.134E  04 
-0.1 26E  02 
-0.134E  04 
-0.137E  02 
-0.134E  04 
-0.1S1E  02 
-0.134E  04 
-0.147E  02 
-0.1 34E  04 
-0.186E  02 
-0.134E  04 
-0.210E  02 
-0.134E  04 
-0.238E  02 
-0.134E  04 
-0.274E  02 
-0.134E  04 
-0.316E  02 
-0.136E  04 
-0.370E  02 
-0.132E  04 
-0.409E  02 
-0.6&9E  03 
-0.957E  00 
-0.734E  03 
-0.U1E  01 
-0.A53E  03 
-0.102E  01 
-0.659E  03 
-0.108E  01 
-0.6S9E  03 
-0.113E  01 
-0.6S9E  03 
-0.119E  01 


82 

26.0 

0.301E  05 

8.8 

0.223E-01 

83 

30.0 

0.28SE  05 

8.8 

0.129E-01 

84 

34.0 

0.268E  05 

8.8 

0.560E-01 

85 

38.0 

0.252E  05 

8.8 

-0. 170E-01 

86 

42.0 

0.236E  05 

8.8 

-0.346E-01 

87 

46.0 

0.220E  05 

8.8 

0.149E-01 

88 

50.0 

0.203E  05 

8.8 

-0.193E-01 

89 

54.0 

0.187E  05 

8.8 

-0.121E  00 

90 

58.0 

0.171E  05 

8.8 

-0.199E-01 

91 

62.0 

0.155E  05 

8.8 

0 . 863E-01 

92 

66.0 

0.13BE  05 

8.8 

-0.463E-02 

93 

70.0 

0.122E  05 

8.8 

-0.450E-01 

94 

74.0 

0.106E  05 

8.8 

-0.239E-01 

95 

78.0 

0.896E  04 

8.8 

-0.360E-01 

96 

82.0 

0.733E  04 

8.8 

-0.305E-01 

97 

86.0 

0.570E  04 

8.8 

-0.105E  00 

98 

90.0 

0.407E  04 

8.8 

-0.953E-01 

99 

94.0 

0.244E  04 

8.8 

0.2S0E  02 

100 

98.0 

0.814E  03 

8.8 

-0.172E  03 

0.301E 

05 

-0.659E 

03 

-0.144E 

02 

-0.126E 

01 

0.285E 

05 

-0.659E 

03 

-0.152E 

02 

-0.133E 

01 

0.269E 

05 

-0.659E 

03 

-0.161E 

02 

-0.141E 

01 

0.2S2E 

05 

-0.659E 

03 

-0.172E 

02 

-0.150E 

01 

0.236E 

05 

-0.659E 

03 

-0.184E 

02 

-0.160E 

01 

0.220E 

05 

-0.659E 

03 

-0.197E 

02 

-0. 172E 

01 

0.204E 

05 

-0.659E 

03 

-0.214E 

02 

-0.186E 

01 

0.187E 

05 

-0.659E 

03 

-0.233E 

02 

-0.202E 

01 

0.171E 

05 

-0.659E 

03 

-0.254E 

02 

-0.221E 

01 

0.1SSE 

05 

-0.6S9E 

03 

-0.280E 

02 

-0.244E 

01 

0.139E 

05 

-0.660E 

03 

-0.314E 

02 

-0.273E 

01 

0.122E 

05 

-0.660E 

03 

-0.3S6E 

02 

-0.309E 

01 

0. 106E 

05 

-0.659E 

03 

-0.410E 

02 

-0.3S6E 

01 

0.900E 

04 

-0.660E 

03 

-0.484E 

02 

-0.420E 

01 

0.739E 

04 

-0.660E 

03 

-0.590E 

02 

-0.511E 

01 

0.S78E 

04 

-0.660E 

03 

-0.756E 

02 

-0.6S3E 

01 

0.418E 

04 

-0.661E 

03 

-0.105E 

03 

-0.901E 

01 

0.260E 

04 

-0.643E 

03 

-0.135E 

03 

-0.140E 

02 

0.U6E 

04 

-0.677E 

03 

-0.516E 

03 

-0.270E 

02 
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Listing  of  Program  on  WES  Computer  for  Generating 
the  Final  Stress  Data  File 


1000C 

1010C  PROGRAM  TO  PRODUCE  STRESS  DATA 

1020C 

1030  CHARACTER  NAME*6»  FNAME*8 

1040  LN  =  10000 

1050C 

L060C  ATTACH  INPUT  STRESS  DATA  FILE 

1070C 

1080  PR  I  NT  »  'ENTER  INPUT  FILENAME* 

1090  READ  100.  NAME 

1100  100  FORMAT  (A6) 

1110  ENCODE  ( FNAME .110)  '/'»  NAME  » 

1120  110  FORMAT  (Al.  A 6.  Al) 

1130  CALL  ATTACH(01.  FNAME  »  3.  0,  1ST) 

1 140C 

1150C  ATTACH  OUTPUT  STRESS  FILE 

1160C 

1170  PRINT  »  'ENTER  OUTPUT  FILENAME* 

1180  READ  100.  NAME 

1190  ENCODE  (FNAME » 1 10 )  '/'.  NAME. 

1200  CALL  ATTACH ( 02 »  FNAME.  3.  0.  1ST) 

1210C 

1220C  ITY=TYPE( ITY=0  FOR  STRESSES  AT  CENTROID) 

1230C 

1240  ITY  =  0 

1250C 

1260C  ITC  =  1  (POSITIVE  STRESSES  ARE  TENSION) 

1270C 

1280  ITC  =  1 

1290C 

1300  WRITE  (2.120)  LN.  ITY.  ITC 

1310  120  FORMAT  (315) 

1320  LN  =  LN  +  10 

1330C 

1340  PRINT.  'ENTER  NUMBER  OF  ELEMENTS' 

1350  READ  130.  NEL 

1360C 

1370  DO  ISO  I  »  1.  NEL 

1380  READ  (1.130)  lELr  XC.  XSTR.  XDUM.  XYSTR 

1390  READ  (1.130)  YC.  YSTR 

1400  130  FORMAT  (V) 

1410  WRITE  (2.140)  LN.  XSTR.  YSTR.  XYSTR 

1420  140  FORMAT  (15.  3E12.3) 

1430  150  LN  -  LN  +  10 

1440C 

1450  STOP 

1460  END 


* 
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LIST  FSTRF1 


10000 

0  1 

10010 

-0.399E 

05 

0.204E  03 

-0.649E 

03 

10020 

-0.382E 

05 

-0.331E  02 

-0.73AE 

03 

10030 

-0.344E 

05 

-0.115E  02 

-0.A53E 

03 

10040 

-0.350E 

05 

0.317E  01 

-0.6S9E 

03 

10050 

-0.333E 

05 

-0.434E  00 

-0.659E 

03 

10040 

-0.317E 

05 

0.605E-01 

-0.459E 

03 

10070 

-0.301E 

05 

0.409E-02 

-0.659E 

03 

10080 

-0.285E 

05 

-0.39AE-02 

-0.659E 

03 

10090 

-0.248E 

05 

-0.136E-01 

-0.659E 

03 

10100 

-0.252E 

05 

-0.315E-01 

-0.439E 

03 

10110 

-0.234E 

05 

0.105E-01 

-0.A59E 

03 

10120 

-0.220E 

05 

0.110E  00 

-0.4S9E 

03 

10130 

-0.203E 

05 

0.453E-01 

-0.659E 

03 

10140 

-0.187E 

05 

-0 . 789E-01 

-0.A59E 

03 

10150 

-0.171E 

05 

-0.755E-01 

-0.659E 

03 

10140 

-0.155E 

05 

0.1S8E  00 

-0.640E 

03 

10170 

-0.138E 

05 

0.131E  00 

-0.659E 

03 

10180 

-0.122E 

05 

-0.170E  00 

-0.660E 

03 

10190 

-0.104E 

05 

0.134E-01 

-0.640E 

03 

10200 

-0.894E 

04 

-0. 459E-01 

-0.660E 

03 

10210 

-0.733E 

04 

-0.192E  00 

-0 . 6&0E 

03 

10220 

-0.570E 

04 

0 . 988E-01 

-0.660E 

03 

10230 

-0.407E 

04 

0.383E  00 

-0.661E 

03 

10240 

-0.244E 

04 

-0.249E  02 

-0.643E 

03 

102S0 

-0.814E 

03 

0.172E  03 

-0.677E 

03 

10240 

-0.130E 

OS 

0.974E  03 

-0.133E 

04 

10270 

-0.127E 

05 

-0.170E  03 

-0.12AE 

04 

10280 

-0.1 22E 

05 

0.103E  02 

-0.134E 

04 

10290 

-0.117E 

05 

-0.417E  00 

-0.134E 

04 

10300 

-0.111E 

05 

0.216E  00 

-0. 134E 

04 

10310 

-0.104E 

05 

-0 . 443E-01 

-0.134E 

04 

10320 

-0.100E 

05 

0.321E-01 

-0.134E 

04 

10330 

-0.949E 

04 

0.130E-01 

-0.134E 

04 

10340 

-0.895E 

04 

-0.212E-01 

-0.134E 

04 

10350 

-0.841E 

04 

-0 . 45BE-01 

-0.134E 

04 

10340 

-0.784E 

04 

-0. 186E-01 

-0.134E 

04 

10370 

-0.732E 

04 

-0. 121E-03 

-0.134E 

04 

10380 

-0.478E 

04 

-0.133E  00 

-0.134E 

04 

10390 

-0.424E 

04 

-0.181E  00 

-0.134E 

04 

10400 

-0.570E 

04 

-0.107E  00 

-0.134E 

04 

10410 

-0.315E 

04 

-0.743E-01 

-0.134E 

04 

10420 

-0.441E 

04 

-0.134E  00 

-0.134E 

04 

10430 

-0.407E 

04 

-0  *  364E  00 

-0.134E 

04 

10440 

-0.353E 

04 

-0.803E-01 

-0.134E 

04 

10450 

-0.299E 

04 

0.464E-01 

-0.134E 

04 

10440 

-0.244E 

04 

-0.722E-01 

-0.134E 

04 

10470 

-0.190E 

04 

0.431E  00 

-0.134E 

04 

10480 

-0.134E 

04 

-0.104E  01 

-0.134E 

04 

10490 

-0.803E 

03 

-0.1 95E  02 

-0.136E 

04 

10500 

-0.24SE 

03 

0.140E  03 

-0.132E 

04 

10510 

0.130E 

05 

-0.976E  03 

-0.133E 

04 

10520 

0.127E 

05 

0.170E  03 

-0.12AE 

04 

10530 

0.122E 

05 

-0.103E  02 

-0. 134E 

04 

10540 

0.U7E 

05 

0.421E  00 

-0.134E 

04 

10550 

0 . 11 IE 

05 

-0.217E  00 

-0.134E 

04 

10540 

0.104E 

05 

0.435E-01 

-0.134E 

04 

10570 

0.100E 

05 

0. 173E-01 

-0.134E 

04 

10580 

0.949E 

04 

0.502E-01 

-0.134E 

04 

10590 

0.895E 

04 

-0.421E-01 

-0.134E 

04 

10400 

0.841E 

04 

-0.991E-01 

-0.134E 

04 

10410 

0.7B4E 

04 

0. 200E-01 

-0.134E 

04 

10420 

0.732E 

04 

0.994E-01 

-0.134E 

04 

10430 

0.478E 

04 

-0.585E-01 

-0.134E 

04 

10440 

0.424E 

04 

-0.110E  00 

-0.134E 

04 

10450 

0.570E 

04 

-0.289E-02 

-0.134E 

04 

10440 

0.515E 

04 

0.406E-02 

-0.134E 

04 

10470 

0.441E 

04 

-0.202E  00 

-0.134E 

04 

10480 

0.407E 

04 

-0.240E  00 

-0.134E 

04 

10490 

0.353E 

04 

-0.582E-01 

-0.134E 

04 

10700 

0.299E 

04 

-0.593E-01 

-0.134E 

04 

10710 

0.244E 

04 

0.428E-01 

-0.134E 

04 

10720 

0.190E 

04 

-0.411E  00 

-0.134E 

04 

10730 

0.134E 

04 

0.980E  00 

-0.134E 

04 

10740 

0.803E 

03 

0.194E  02 

-0.136E 

04 

10750 

0.245E 

03 

-0.140E  03 

-0. 132E 

04 

10700 

0.399E 

05 

-0.206E  03 

-0.4A9E 

03 

10770 

0.382E 

05 

0.331E  02 

-0.736E 

03 

10780 

0.344E 

05 

O.USE  02 

-0.6S3E 

03 

A23 


10790 

0.330E 

03 

-0.317E  01 

-0.659E 

03 

loeoo 

0.333E 

03 

0.434E  00 

-0.659E 

03 

10810 

0.317E 

OS 

-0.472E-01 

-0.659E 

03 

10820 

0.301E 

03 

0.223E-01 

-0.639E 

03 

10830 

0.285E 

03 

0.129E-01 

-0.659E 

03 

10840 

0.268E 

05 

0.560E-01 

-0.6S9E 

03 

10830 

0.232E 

03 

-0.170E-01 

-0.639E 

03 

10860 

0.236E 

03 

-0.346E-01 

-0.6S9E 

03 

10870 

0.220E 

03 

0.149E-01 

-0.6S9E 

03 

10880 

0.203E 

03 

-0.193E-01 

-0.6S9E 

03 

10890 

0.187E 

03 

-0.121E  00 

-0.659E 

03 

10900 

0.171E 

03 

-0.199E-01 

-0.639E 

03 

10910 

0.155E 

05 

0.863E-01 

-0.6S9E 

03 

10920 

0.138E 

05 

-0.463E-02 

-0.660E 

03 

10930 

0.122E 

03 

-0.450E-01 

-0.660E 

03 

10940 

0.106E 

05 

-0.239E-01 

-0.639E 

03 

10930 

0.896E 

04 

-0.360E-01 

-0.660E 

03 

10960 

0.733E 

04 

-0.305E-01 

-0.660E 

03 

10970 

0.570E 

04 

-0.103E  00 

-0.660E 

03 

10980 

0.407E 

04 

-0.953E-01 

-0.661E 

03 

10990 

0.244E 

04 

0.230E  02 

-0.643E 

03 

11000  0.814E  03  -0.172E  03  -0.677E  03 


1.  The  following  problem  is  for  the  cross  section  of  a  lock  wall.  This 
same  problem  has  been  used  in  the  WES  report  (in  preparation)  "Case  Study  of 
Six  Major  General-Purpose  Finite  Element  Programs."  The  dimensions  of  the 
problem  as  well  as  contour  plots  of  the  stresses  are  given  in  that  report. 

The  cross  section  has  only  one  material  type.  The  purpose  of  the  problem  is 
to  demonstrate  the  sequence  of  commands  that  might  be  used  for  a  problem  of 
this  type. 

2.  Data  files  GEO  and  STRETM  were  created  in  the  same  manner  as  those 
described  in  Appendix  A.  The  file  GEO  contains  the  necessary  geometry  (node 
and  element)  data  and  the  file  STRETM  contains  the  stress  data. 

Defining  files 

ENTER  GEOMETRY  DATA  FILE  NAME ' 

-GEO 

ENTER  STRESS  DATA  FILE  NAME: 

“STRETM 

ENTER  OUTPUT  DATA  FILE  NAME* 

-OUTEX 

Command  for  plotting  grid 


COMMAND? 
=P  G 


Defining  sections 


couhanp* 

•  s  1 

COWWND? 

•  581 
COMMAND? 

•  S  3  S  18  8?  ta 

command? 

•  S-4  38  18  11«  18  1 
COMMAND? 


Section  1  (cross  hairs  used  for  section  coordinates) 


Section  2: 


Section  3: 


Section  4: 


Consider  only  material  type  //I  (cross  hairs 
used  for  section  coordinates) 

XI  =  18.  Y 1  =  12.  X2  =  27. 

Y2  =  12. 

XI  =  38.  Y1  =  12.  X2  =  50. 

Y2  =  12.  Material  Type  =  1 


Output  for  section  4 


!«:  v|>  •  (twill*.) 

"""■S3S :  :&&*•’ 

MBIT  •  -  .Mim 
THRUST  -  -.SISK'S 


.It.) 


KCT10H  NS.  4 


Output  file 

3.  The  following  is  the  output  file  for  this  problem.  It  contains 
two  lines  of  data  for  each  section  which  has  been  displayed.  The  section 
number,  XI,  Yl,  X2,  Y2,  and  location  of  the  neutral  axis  are  given  on  the 
first  line.  The  second  line  contains  values  for  shear,  moment,  and  thrust 


•LIST  OUTEX 

1  33.87 
2.1S992H:  26 

2  -2.74 
2.139242C  M 

3  I.N 
2.293192C  22 

4  32.22  _ 

MSS24K  M 


-I.lt  34. M 

-2.2S9741C  K 
2.13  112.14 

4. 34 TBS 1C  27 
12. M  27.29 

4.173724C  M 

12.22  112.22 

-2. 8921 ISC  22 


12.22 

-2.224I9SC 

12.22 

-2.I1224SC 


34.22 

2.27 

2S 

92.22 

22 

2.32 

19.99 

12.22 

22 

92.92 

12.22 

22 
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